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AWARDS ABSTRACT 


This invention relates to a method for digital reduction of synthetic 
aperture multipolarized radar data while still allowing full polarimetric 
utility of the data. The technique results in about 12.8-fold reduction in 
volume to allow multipolarized radar data to be used by a university and small 
industry users. 

FIG. 1 illustrates a flow chart of a prior art process for generating a 
synthetic polarization image from a set of imaging radar polarimeter data, FIG. 
2 illustrates a flow chart of a first new method (scattering matrix approach) 
for generating a synthetic polarization image with data volume reduction, and 
FIG. 3 illustrates a flow chart of a second new method (Stokes matrix approach) 
for generating a synthetic polarization image with data reduction. These flow 
charts define the alternative methods of the invention. In the first embodi- 
ment of FIG. 2, the scattering matrices of four consecutive along-track pixels 
are "averaged," while in the second embodiment of FIG. 2 the Stokes matrices of 
four consecutive along-track pixels are "averaged." In both cases, the average 
matrices are stored as a reduced data set. Each approach in terms of data 
volume reduction and in terms of errors introduced in the synthesized images 
substantially reduces the size of memory required for each of the images, and 
the number of operations is also reduced, thus reducing data processing time. 
This is accomplished in the approach of FIG. 2 by compressing scattering matri- 
ces, and in the approach of FIG. 3 by compressing Stokes matrices. The reduced 
data set is smaller in the first approach than that created by compressing 
Stokes matrices. However, greater error is introduced by compressing scattering 
matrices: these errors may be as great as 10$ to be compared to typically less 
than 10 - 3 when the compression algorithm operates on the Stokes matrices. FIG. 
4 is a typical histogram of the Stokes matrix elements F mn used as an example 
in determining the preferred method for optimizing the Stokes matrix method of 
data reduction illustrated in FIG. 3- 

The invention resides in the "averaging" technique of the two alterna- 
tive methods to reduce data corresponding to four look pixels. 
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DATA VOLUME REDUCTION FOR 
IMAGING RADAR POLARIMETRY 


Ori g in of the Inv ention 

5 The invention described herein was made in the 

performance of work under a NASA contract, and is 
subject to the provisions of Public Law 96-517 (35 
USC 202) in which the Contractor has elected not to 
retaintitle. 

10 

Tec hnical Fie ld 

This invention relates to a method for digi- 
tal reduction of synthetic aperture multipolarized 
radar data while still allowing full polarimetric 
15 utility of the data. The technique results in about 
12.8-fold reduction in volume to allow multipolarized 
radar data to be used by a university and small 
industry users. 


20 Bac kgrou nd Art 

Recent radar measurements show that different 
ground locations respond differently when the polari- 
zation of either the receiving or the transmitting 
antenna is varied (see Dino Guili, "Polarization 
25 diversity in radars," Proceedings of the IEEE, Vol. 
74, No. 2, February 1986). In 1985, a polarimeter 
was flown on a NASA aircraft, recording data which 
were subsequently processed and stored at the Jet 
Propulsion Laboratory. One difficulty inherent in 
30 this experiment was the large storage necessary for 
each of the images; only a limited number of data 
sets could be stored. Furthermore, the computer time 
necessary to synthesize a picture using an arbitrary 
transmit and receive polarization is mostly devoted 
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to cumbersome data transfers. Clearly, the volume of 
data needs to be reduced in order to provide the user 
with a more flexible investigation tool, yet the data 
volume reduction must not impact the noise level by 
introducing additional error. Some theoretical con- 
cepts will first be presented for clarity and review, 
as well as a description of the operations necessary 
to generate a synthetic polarization image from the 
original data sets. Then two new methods of data 
compression will be described. 

State m ent of th e Invention 

In a first embodiment of the invention, the 
scattering matrices of four consecutive along-track 
pixels are "averaged,'* while in a second embodiment, 
the Stokes matrices of four consecutive along-track 
pixels are "averaged." In both cases, the average 
matrices are stored as a reduced data set. Each ap- 
proach in terms of data volume reduction and in terms 
of errors introduced in the synthesized images sub- 
stantially reduces the size of memory required for 
each of the images, and the number of operations is 
also reduced, thus reducing data processing time. 
This is accomplished in one approach by compressing 
scattering matrices, and in the other by compressing 
Stokes matrices. The reduced data set is smaller in 
the first approach than that created by compressing 
Stokes matrices. However, greater error is introduced 
by compressing scattering matrices: these errors may 
be as great as 10$ to be compared to typically less 
than 10~3 when the compression algorithm operates on 
the Stokes matrices. 
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B ri ef Desc ription o f th e D ra wi ngs 

FIG. 1 illustrates a flow chart of a prior art 
process for generating a synthetic polarization image 
from a set of imaging radar polarimeter data. 

5 FIG. 2 illustrates a flow chart of a first new 

method (scattering matrix approach) for generating a 
synthetic polarization image with data volume reduc- 
tion. 

FIG. 3 illustrates a flow chart of a second 
10 new method (Stokes matrix approach) for generating a 
synthetic polarization image with data reduction. 

FIG. is a typical histogram of the Stokes 
matrix elements F mn used as an example in determining 
the preferred method for optimizing the Stokes matrix 
15 method of data reduction illustrated in FIG. 3- 


20 


25 


Detail ed D escrip ti on of the Inven tion 

The expression for the electric field of an 
electromagnetic wave propagating along the z axis is: 


E ( z , t ) = Re 


Er ( z , t ) 
Ey ( Z | t ) 


Re ( h e i(kz - wt) ) 


(1 ) 


where, h 


a v e 


-15 


exp(i6 H ), and 6 = 6r-S v 


where h is the polarization vector, Er and E v are the 
horizontally and vertically polarized components of 
30 the electric field with respective magnitudes a H and 
a v and arguments 6 r and 5y. 

Another way to describe the same electric 
field is given by the Stokes parameters, defined as 
foil ows : 
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a H " a V 
2a H a v cos ( 6 ) 

2a H a v si n( 5 ) 


( 2 ) 


The Stokes vector is the vector whose components are 
the four Stokes parameters. If Gq 2 = G .| 2 + G 2 2 + 

G^ , the electromagnetic wave is said to be fully 
polarized . 

If the backscatter i ng phenomena is assumed to 
be linear, isotrospic and homogeneous, the back- 
scattered wave polarization vector h s may be ex- 
pressed as 


h s = Sh fc where S 


S S 

xx xy 

c c 

yx yy 


(3) 


where h t is the polarization vector of the incident 
wave or transmitting antenna and S, the scattering 
matrix, is a 2x2 complex matrix. If h p is the po- 
larization vector of the receiving antenna, the 
complex ampl i tude' - of the received signal is: 


V = h p T Sh t (4) 

Thus, once the scattering matrix is known, a 
synthesized response may be computed for any desired 
configuration of antenna polarization states defined 
by h p and h fc . 

Similarly, the 4x4 real Stokes matrix F 
relates the power of the signal received to the 
Stokes vectors defining the polarization state of the 
receiving and transmitting antennas, G_ and Gj. . 
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If the reciprocity principle applies, then the 
scattering matrix is symmetrical and so is the Stokes 
matrix. In that case, the Stokes matrix can be 
derived from the scattering matrix as follows (the 
general relation can be found in "Light scattering by 
small particules," H.C. Van de Hulst, Dover publica- 
tion) : 


xxxx 

- q c * 

S XX S XX 

^xxxy 

" S xx S xy* 

yyyy 

S yy S yy* 

J xyyy 

= S xy s yy * 

xy xy 

= S xy S xy* 

J xxyy 

" S xx S yy* 


where the asterisk indicates the conjugate of the 
signal indicated. 


F H " °' 25 ^ J xxxx +2 J xy xy + **yyyy ^ 

F 12 = °* 25 (J xxxx- J yyyy) 

F 1 3 = 0,5 Re ^ J xxxy +,J xyyy ^ 

F 14 " °- 5 Im ^ -J xxxy~ J xyyy ^ 

F 44 = 0.5 Re (-J xxyy +J xyxy > 

There is enough information in either the 
Stokes matrix, F, or the scattering matrix, S, to 
synthesize any polarization configuration. The 

scattering matrix yields a synthesized signal whose 
power is identical to the power obtained from the 
Stokes matrix as long as the waves are fully polar- 
ized. 

A wave is fully polarized if it can be ex- 
pressed as the superposition of a horizontally polar- 
ized part (HP) and a vertically polarized part (VP). 
The polarization vector describes the electric field 
as a combination of a vertical part and a horizontal 
part. Therefore this representation is always refer- 


f 23 

* 0.5 

Re 

( J xxxy“ J xyyy ^ 

F 24 

= 0.5 

Im 

^ -J xxxy +J xyyy ^ 

F 33 

= 0.5 

Re 

^ J xxyy +J xyxy ^ 

f 34 

= 0.5 

Im 

^~ J xxyy ^ 


( 6 ) 


L 


J 
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ring to a fully polarized wave. Since the scattering 
matrix is the linear operator associated with a 
polarization vector, the scattering matrix represen- 
tation of a scatterer assumes that this scatterer 
5 cannot introduce any diffuse component (part of a 
wave which is not polarized) in the backscattered 

wave. The backscattered wave has to be fully polar- 

ized because it is written as a polarization vector. 
The Stokes vector representation, as opposed to the 

10 polarization vector representation, allows for a 
diffuse component which can be estimated by the 
difference g q 2- ( G i 2+G 2 2+G 3 2 ^ (see H.C. Van de Hulst, 
supra). Consequently, the Stokes matrix representa- 
tion of a scattering surface can include this depol- 
15 arization phenomena. 

Three methods will now be described with a 
complete set of data from an imaging radar polari- 
meter flown in 1985 on a NASA CV990 aircraft. Out of 
the subsequent data set, three images were chosen 
20 which feature a variety of targets for experimental 

testing of an original (prior art) process illus- 

trated in FIG. 1 and two new processes, one using a 
scattering matrix approach illustrated in FIG. 2 and 
the other using a Stokes matrix approach illustrated 
25 in FIG. 3. The original process requires about 128 

megabytes of synthetic aperture image data storage 
for one scene (frame). The complete data set is 

comprised of 4.2 million one-look pixels, 1024x4 
points in the along-track direction and 1024 pixels 
30 in the range direction. For each pixel, four complex 
elements are stored per scattering matrix. A, B, C, 
and D in blocks 10, 11, 12, and 13 represent four 

consecutive along-track pixels. Since a complex 

number is 8 bytes long, the storage requirements are 
35 128 megabytes. If the reciprocity principle is 
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assumed, then the scattering matrix should be symmet- 
rical and only the three different elements of the 
scattering matrix need be stored. In fact, this 
property is assumed when calibrating the data. 

5 In the original process, it was necessary to 

generate one synthetic polarization image from this 
set of data in the following manner. The polariza- 
tion vectors of both the transmitting and the receiv- 
ing antenna (not shown) of the experiment are 
10 chosen. The received signal is synthesized for each 
pixel using equation ( 4 ) , as indicated by blocks 15 , 
16 , 17 and 18 . The resulting powers are then com- 
puted in blocks 19 , 20 , 21 and 22 . The last step 

consists of averaging the computed power of 4 con- 
15 secutive along-track points in block 23 . This opera- 
tion reduces statistical variations and makes the 
element of resolution approximately square. Each 
resulting averaged pixel now has an associated power 
and the image is ready for display, as indicated in 
20 block 24, namely as intensity relative to one four- 
look pixel . 

This process of synthesizing an image of arbi- 
trary polarization requires about twenty 
minutes on a VAX 785 computer. Future users 
25 of the data will want to interact quickly with the 

picture making a faster process highly desirable. 
Reducing the size of the data set decreases the 
processing time. The consequently reduced storage 
requirements also allow the storage of more data sets 
30 where only a few were previously possible. 

Two different data reduction methods have been 
developed for reduced storage requirements using 
different approaches referred to hereinbefore as the 
scattering matrix approach and the Stokes matrix 
35 approach. In each case, a matrix corresponding to 
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four along-track consecutive pixels A, B, C and D is 
stored as the result of an "averaging" operation done 
directly on the matrices, and not, as in the original 
process just described above, done on the synthesized 
5 signal power. The two new methods are detailed in the 
following sections, and a comparative study of the 
results is presented in a subsequent section. In 
each approach, the four complex elements A, B, C and 
D of the scattering matrix are stored as indicated at 
10 the top of FIGs. 2 and 3 in blocks 10, 11, 12 and 13* 

The Scattering Matrix Approach 

An "average" of the scattering matrices is 
computed rather than an average of the synthesized 
15 signal powers: the "average" scattering matrix corre- 

sponding to one four look pixel is then the only 
information needed to synthesize a picture. FIG. 2 
is a flow chart which describes the process. For 
convenience in understanding this first embodiment 
20 and the second to be described afterwards, the same 

reference numerals are employed for the same elements 
or blocks as used in FIG. 1. First an average scat- 
tering matrix is computed in block 25 for every set 
of four consecutive points in the azimuth direction. 
25 This matrix is stored in 8 bytes of computer memory 
represented by block 26. Eight megabytes of memory 
are then required to store the entire reduced data 
set for one image . 

Each time it is desired to generate a picture 
30 from a reduced data set 27, the first step is to 

choose the transmit and received polarization vectors 
in block 14. The synthesized signal is then computed 
using Equation (4) and its power from P = VV*, but in 
this case the reduced data set is of averaged scat- 
35 tering matrices of four-look pixels, as indicated by 
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blocks 28, 29 and 30, thus reducing the number of 

operations required to synthesize a complete image. 
The process is faster than before because less, data 
is being manipulated and fewer operations are being 
5 performed. Data reduction can only be effectively 

achieved by carefully choosing the "averaging" me- 
thod. The scattering method will now be described in 
more detail. 

Let A, B, C and D be the four scattering 
10 matrices to be "averaged," corresponding to four 
consecutive points in the azimuth direction and let S 
be the resulting matrix. 


A = 

A xx 

A xy 

B = 

B xx 

B xy 


A yx 

Ayy 


B y x 

B yy 

C = 

C X x 

C xy 

D = 

c xx 

o 

X 


C yx 

c yy 


C yx 

c yy 


In the output image, each pixel is associated 
simply with an intensity, and the absolute phase 
information of the four received signals can be 
discarded, although the relative phases are yet re- 
25 quired. Therefore, the absolute phase of the scat- 
tering matrix can be set arbitrarily. The power in 
each element of the "averaged" matrix is computed as 
follows : 

30 l S xxl = s ^ r ^ A xx A xx* + B xx B xx* + t 'xx^xx* + D xx D xx*^ 

| S X y | = S 9 r t ( A X y A x y * + B xy B xy *+ C x y C x y *+ D x yD x y * ) 

|Sy X | = Sqr t ( Ay x Ay x * + By x By x *+ Cy X Cy x *+ Dy X Dy x * ) 

| S yy | = Sqrt(AyyAyy* + ByyB y y* + CyyCyy* + DyyDyy*) 


( 7 ) 
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where "sqrt" signifies "the square root of In 

the following, the phase of S xx is set to zero. As 
discussed before, this is legitimate as long as the 
phases of the other terms relative to S xx are pre- 
5 served. 

This method for calculating the phase puts 
more weight on the vectors of large amplitude which 
are presumably less sensitive to noise. The informa- 
tion to be coded consists of four amplitudes and 
10 three phases and can be stored in 8 bytes as follows: 

2 bytes for the maximum amplitude stored as a 
mantissa and an exponent. 

15 3 bytes for the three remaining amplitudes. 

3 bytes for the three phases. 

After reduction, the data consists of 
20 1 024x1 02 M pixels and 8 bytes per pixel, resulting in 

an overall reduction ratio of 16. The reduced data 
set may then be processed to synthesize a complete 
image, using one "averaged" scattering matrix per 
pixel, as indicated in block 28 and then proceeding 
25 in blocks 29 and 30 by computing only one signal 

using Equation (4) and then computing power by the 
equation P = VV* to produce the same signal in block 24 
as in the original method. Thus, in reducing data by 
averaging the four scattering matrices in block 25 , 
30 there is only one signal processing path in synthe- 

sizing a complete image, instead of the four in the 
original method illustrated in FIG. 1. 
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The Stokes Matrix Approach 

FIG. 3 is a flow chart which describes the 
second new method for data volume reduction. Using 
Equation (5), each scattering matrix A, B, C and D in 
blocks 10 , 11 , 12 and 13 is transformed into its 
corresponding Stokes matrix F^, F 2 , F^ and F ^ in 
blocks 31. 32, 33 and 3$. One interesting property 
associated with the Stokes matrix is that it yields 
directly the power given the Stokes vectors of both 
the receiving and the transmitting antennas (Equation 
7). Therefore the two following processes are equi- 
valent; adding the synthesized power of the signal 
scattered from 4 different areas or adding the 4 
Stokes matrices characteristic of the four areas, and 
then computing the power of the resulting signal from 
this composite matrix. Let Fj, F 2 , F^, F 14 be the 
Stokes matrices associated with four consecutive 
pixels and, let G t and G p be the Stokes vectors of 
the transmitting and receiving antennas, respective- 
ly. The powers corresponding to the phase matrix 
Fj are given by: 


p l - °r TF 

i G t 


( 8 ) 

The four-look averaging process 

consists 

of 

adding the powers of four 

consecutive 

pixels 

as 

indicated in block 35. The 

pressed as follows: 

resulting P 

can be 

ex- 

P - Z P A - Z G r T F i G t 

= G r T (EF i )G t 


(9) 


Therefore it is possible to add the four 
Stokes matrices corresponding to four consecutive 
points in the azimuth direction to form a single 
Stokes matrix, as indicated in FIG. 3 by block 36 . 
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This process is equivalent to a four-look averaging 
operation. The resulting Stokes matrix correspond- 
ing to a four-look pixel is then stored in block 37 
in a compressed form as the reduced data set. 

5 To form an image from this reduced data set 

stored in block 37, the transmit and received Stokes 
vectors are computed (blocks 38) from the correspond- 
ing polarization vectors I 1 ! in Equation (2). For 
each pixel, the Stokes matrix (block 39) is selected 
10 from the reduced data set 37 and the synthesized 

signal power is obtained in block 40 through Equation 
( 7 ) . 

Each Stokes matrix in block 36 is a 4x4 sym- 
metrical real matrix. It consists of 10 distinct 
15 elements, nine of which are independent (Equation 5). 
The tenth, F 22 » can be computed from three of the 
nine elements as follows: 

F 22 F 1 I ” F 33 F 44 (10) 

20 

F-j i is the largest element and is always 
positive. It is coded in two bytes as a mantissa and 
an exponent. The eight other independent elements 

are scaled by F 11 and coded in 1 byte each. There- 
25 fore the complete phase matrix can be stored using 10 
bytes. The complete data set for one image 

(1024x1024 pixels) is now stored in 10x1024x1024 
bytes. The overall reduction factor from the origi- 
nal data is 12.8. 

30 

Comparative Study of the Results ' 

of the Different Algorithms 

A testing procedure necessary to validate the 
different algorithms and to compare their perform- 
35 ances compares a polarization signature obtained from 
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the original data set with the corresponding polari- 
zation signature computed from the reduced data set. 
A polarization signature is a surface in which height 
corresponds to the intensity of the signal received 
5 by an antenna transmitting a polarization defined by 
the x and y coordinates of the point on the surface. 
The output of this program is: 


10 


15 


20 


25 


Error 


V //(I-I) 2 slnxdxd4) 

m 2 sl \V* 


(ID 


where I is the intensity of the pixel computed from 
the original data and I is the intensity computed 
from the reduced data set. The integration is com- 
puted over the set of all possible polarizations for 
the emitting antenna assuming the receiving antenna 
is the same as the emitting antenna (copolarized 
error) or assuming that the receiving and emitting 
antennas are orthogonal or crosspolari zed (cross- 
polarized error ) . 

The test program outputs shown below indicate 
the difference between the original polarization 
surface and the one computed from the reduced data 
set for both reduction methods. 


ERROR (Eq. 11) 

Forest area 

Urban area 

Ocean area 

copol 

signa 

Scattering 

matrix 

method 

7.19 E-2 

4.10 E-2 

8.80 E-2 

Stokes 

matrix 

method 

4.06 E-4 

3.40 E-4 

3.14 E-4 
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ERROR (Eq. 11) 

Forest area 

Urban area 

Ocean area 

cross- 

pol 

signa 

Scattering 

matrix 

method 

1.20 E-1 

7.20 E-2 

1.20 E-1 

Stokes 

matrix 

method 

4.18 E-4 

3.09 E-4 

4.25 E-4 


The Stokes matrix method yields smaller error 
10 than the scattering matrix method. The four-look 

averaging process on the original data set corre- 
sponds to adding the power of four consecutive pixels 
and therefore is completely equivalent to adding the 
four corresponding Stokes matrices. The scattering 
15 matrices yield the complex field, not the power of 

the received signal and "averaging" four consecutive 
scattering matrices is not equivalent to a four-look 
average. The Stokes matrix representation allows for 
a diffuse component while the scattering matrix 
20 representation ignores it. For this reason the 

Stokes matrix compression introduces smaller error in 
the data. The compression ratio is better in the 

case of the scattering method but the resulting 
errors are not acceptable. The Stokes matrix algo- 
25 rithm was chosen and can be further improved as 

foil ows . 

Optimization of the Stokes Matrix Method 

The error computed above can be somewhat 
30 reduced by more carefully choosing the storage proce- 

dures. The method preferred is based on a probabil- 
istic knowledge of the relative magnitudes of terms 
in the Stokes matrix. 

As discussed previously, F 1 1 , the largest 
35 element of the Stokes matrix, is coded in two bytes. 
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The first one corresponds to the exponent EXP and the 
second one corresponds to the mantissa MAN. 

EXP = E(ln(F 1 1 ) / ln2) , 

5 ( 12 ) 

MAN = F n /2 EXP 

where E(x) is the largest integer smaller than x. 

Since MAN is' always less than 2 and greater or 
10 equal to 1, better precision can be achieved by 

storing MAN as: 

by te ( 2 ) - E( 254*(MAN-1 . 5) ) (13) 

15 where byte(2) will range between -128 and 127. EXP 

is stored directly as one byte: 


byte(1) = EXP 


(14) 


20 The eight remaining elements are first scaled 

by F n 



F ij /F ll 


(15) 


25 The absolute value of each P mn , where the 

subscripts mn designate a particular Pjj» is always 
less than 1. Two different strategies are possible 
for the storage of each term: storing P mn directly 

or storing its square root. If P mn is digitized in 
30 1 byte, the resulting error is ±2 - ®. If /P mn is 

digitized, the error on /P mn is ±2~®, therefore the 
resulting P mn is: 


(/P +2 - ®) 2 =P ± /P 2“^ 

VT mn~ 1 mn 1 * mn 
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Comparing 2“® and 2~^/P mn shows that taking 
the square root yields better results as long as /P mn 
is less than 0.5, i.e., P mn <0.25. A typical histogram 
of each P mn from actual data collected over the San 
5 Francisco Bay area is presented in FIG. 4. p i ^ » P 33 * 
P ij ij are mostly larger than 0.25 while P 12 » P 13’ P 14' 
p 2 s . p 24 » p 34 tend to be smaller. 

In the case where the square root of P is 

stored, 

10 

byte(i) = E(sign(P mn )*/T p ~nT* 12 7 ) (16) 

where E(x) is the largest integer smaller than x and 
sign (x)=-1 when x is negative and +1 otherwise. 
15 When P is stored directly, 

byte(i) - E(127*P mn ) (17) 

Four different algorithms, each corresponding 
20 to different numbers of square root elements were 


tested and the table below shows the differences 
between them. 



Algo. A 

Algo. B 

Algo. C 

Algo. D 

P 1 2 

normal 

normal 

normal 

normal 

P 1 3 

normal 

normal 

Square rt 

Square rt 

P 1 4 

normal 

Square rt 

Square rt 

Square rt 

P 23 

normal 

Square rt 

Square rt 

Square rt 

P 2 4 

normal 

normal 

Square rt 

Square rt 

P 33 

normal 

normal 

normal 

normal 

p 3 4 

normal 

normal 

normal 

Square rt 

P 4 4 

normal 

normal 

normal 

normal 
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For example, the square root of P 2 4 is stored 
for algorithm C and D according to Equation (16). In 
algorithm A and B, P 2 n is stored directly (Equation 
( 17 )). 

The tests conducted on these algorithms 
consist of comparing the signature surfaces of a 
specific area, generated from the studied algorithm 
to the signature surfaces generated from the 
unreduced set of data. The test results, which are 
error percentages, are presented below. 


error 

(Eq. 11) Algo. A 


Algo. B Algo. C 


Algo. D 


forest 

CO 

4.06 E-4 

4.46 E-4 

2.80 E. 4 

2.58 E-4 

area 

cross 

4.18 E-4 

3.71 E-4 

4.11 E-4 

3.48 E-4 

urban 

CO 

3.40 E-4 

3.95 E-4 

3.23 E-4 

3.25 E-4 

area 

cross 

3.09 E-4 

2.10 E-4 

2.13 E-4 

2.24 E-4 

ocean 

area 

CO 

3.14 E-4 

3.07 E-4 

2.08 E-4 

2.09 E-4 

cross 

4.25 E-4 

4.29 E-4 

2.51 E-4 

2.51 E-4 


The square root operation requires more proc- 
essing when reducing the data as well as when gener- 
ating a picture. A trade off has to be made between 
better precision and a faster method. Algorithm C 
was finally chosen. The errors introduced by algo- 
rithm C are smaller than when algorithm B is used. 
Algorithm D is more complex and its results are 
comparable with those of algorithm C. 

The data reduction algorithm scheme finally 
adopted involves four square roots (Algorithm C): 
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byte (3) = E(127*P 12 ) 
byte (4) = E(127*sign(P 13 )*/JP73T) 
byte (5) - E( 1 27 * s i gn ( P 1 „ ) VJF^T ) 
byte (6) = E( 1 27 * s i gn ( P 2 3 ) VJp^ 3 T ) (18) 

5 byte (7) = E( 1 27 * s i gn ( ? )VJP 2 J) 

byte (8) = E(127*P 33 ) 
byte (9) = E(127*P 34 ) 
byte (10)= E(127*P 35 ) 

10 Reconstruction of the Stokes Matrix 

To reconstruct the Stokes matrix from the 
reduced data, the following operations are required: 

F ( 1 , 1 ) = (byte(2) /254+1 .5)2 byte(1 > 

15 F ( 1 ,2) = byte(3)*F(1 ,1 j/127 

F( 1 , 3) - sign(byte(4))*[byte(4)/l27] 2 *F( 1 , 1) 

F ( 1 , 4 ) = sign(byte(5) ) *[byte(5)/1 27] 2 *F( 1 , 1 ) 

F ( 2 , 3 ) = si gn ( by t e ( 6 ) ) *[byte(6)/1 27 ] 2 *F ( 1 , 1 ) 

F ( 2 , 4 ) = sign(byte(7) ) *[ by t e ( 7 ) / 1 27] 2 *F( 1 , 1 ) 

20 F ( 3 . 3 ) = byte(8)*F( 1 , 1 )/127 

F ( 3 , 4 ) = byte(9) *F( 1 , 1 )/1 27 (19) 

F( 4, 4) = byte(1 0)*F( 1 , 1 ) / 1 27 
F(2, 2) = F ( 1 , 1 ) - F( 3. 3) " F ( 4 , 4 ) 

25 Conclusion 

A method of reducing required storage space 

for mul t i p ol ar imet r i c synthetic aperture radar data 

has been developed which preserves signal integrity. 

The data reduction operation reduces the required 

30 storage space by 12.8 and speeds the image synthesis 

process by a factor of 10 in time. The errors intro- 

-4 

duced in the output images are of the order of 10 
Thus operations with polarimetric data are greatly 
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facilitated and are within reach of even small re- 
search groups. 

Although particular embodiments of the inven- 
tion have been described and illustrated herein, it 
is recognized that modifications and variations may 
readily occur to those skilled in the art. Conse- 
quently, it is intended that the claims be inter- 
preted to cover such modifications and variations. 
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DATA VOLUME REDUCTION FOR 
IMAGING RADAR POLARIMETRY 


AB STR ACT OF THE DISCLOSURE 

Two alternative methods are disclosed for 
digital reduction of synthetic aperture multipolar- 
ized radar data using scattering matrices, or using 
Stokes matrices, of four consecutive along-track 
pixels to produce "averaged" data for generating a 
synthetic polarization image. 


